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Appendix  to  Final  Report: 
Statement  of  the  problem  studied,  and 
summary  of  the  most  important  results 

Our  project  aimed  to  develop  and  apply  computational  methods  to  investigate  corre¬ 
lation  effects  in  atomic  and  molecular  systems,  and  in  ultracold  atomic  gases.  Below  the 
most  important  results  are  summarized.  The  twelve  (12)  journal  publications  listed  ear¬ 
lier  describe  the  findings  in  detail;  the  work  was  also  presented  in  about  25  invited  talks 
at  international  meetings  [including  at  annual  meetings  of  the  American  Physical  Society 
(APS)  and  the  American  Chemical  Society  (ACS)],  workshops,  and  conferences,  and  also  as 
seminars/colloquia  at  academic  instutitions  around  the  world. 

A  Correlation  effects  in  the  ground  state  of  trapped 
atomic  Bose  gases. 

We  studied  the  effects  of  many-body  correlations  in  trapped  ultracold  atomic  Bose  gases  [1] . 
We  calculate  the  ground  state  of  the  gas  using  a  our  ground-state  auxiliary-field  quantum 
Monte  Carlo  (QMC)  method  [2].  We  examine  the  properties  of  the  gas,  such  as  the  energet¬ 
ics,  condensate  fraction,  real-space  density,  and  momentum  distribution,  as  a  function  of  the 
number  of  particles  and  the  scattering  length.  We  find  that  the  mean-field  Gross-Pitaevskii 
(GP)  approach  gives  qualitatively  incorrect  result  of  the  kinetic  energy  as  a  function  of  the 
scattering  length.  We  present  detailed  QMC  data  for  the  various  quantities,  and  discuss  the 
behavior  of  GP,  modified  GP,  and  the  Bogoliubov  method  under  a  local  density  approxima¬ 
tion. 


B  Phaseless  auxiliary-field  method  —  quantum  chem¬ 
istry  applications  and  benchmark. 

Accuracy,  or  predictive  power,  is  a  key  requirement  of  a  successful  paradigm  for  strongly 
correlated  systems.  The  GTO  formulation  and  chemistry  applications  have  been  crucial  in 
establishing  the  phaseless  AFQMC  method  as  a  general  approach  that  can  now  be  applied 
to  correlated  systems.  To  benchmark  the  method,  we  have  to  date  applied  it  to  close  to  100 
systems,  including  first-  and  second-row  molecules  [3,  4],  transition  metal  oxide  molecules  [5], 
simple  solids  [6,  7],  post-d  elements  [8],  and  van  der  Waals  systems  [9].  Figure  1  illustrates 
the  results  on  molecules  using  both  planewave  and  GTOs.  In  these  calculations  we  have 
operated  largely  in  an  automated  mode,  inputting  only  the  DFT  or  HF  solutions.  The 
method  demonstrated  excellent  accuracy,  consistently  able  to  correct  errors  in  the  mean- 
field  trial  wave  function.  In  molecules,  we  find  that  the  accuracy  of  the  phaseless  AFQMC  is 
comparable  to  the  coupled-cluster  method  CCSD(T),  the  gold  standard  in  chemistry  [10,  11], 
near  equilibrium  geometry.  [The  computational  cost  of  CCSD(T)  scales  as  N7  with  system 
size,  while  QMC  is  N3-N4.] 
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Figure  1:  Calculated  binding  energies  [15,  5, 
4,  16]  of  molecules  compared  with  experimental 
values.  The  discrepancy  between  theory  and  ex¬ 
periment  is  plotted.  The  AFQMC  is  fed  a  trial 
wave  function  to  start,  which  is  taken  directly 
from  DFT  [with  either  local-density  approxima¬ 
tion  (LDA)  or  generalized-gradient  approxima¬ 
tion  (GGA)  functionals]  or  Flartree-Fock  ( H F) . 
The  corresponding  DFT  or  HF  results  are  also 
shown.  As  can  be  readily  observed,  the  AFQMC 
results  are  in  excellent  agreement  with  experi¬ 
ment  and  significantly  improve  upon  the  values 
from  DFT  and  H F. 


C  Electron  correlation  effects  from  bond  stretching  and 
bond  breaking. 

Bond  stretching  in  molecules  and  solids  mimics  different  levels  of  particle  correlation  and 
provides  a  challenging  testbed  for  approximate  many-body  computational  methods.  With 
single  Slater  determinants  from  unrestricted  Hartree-Fock  (UHF)  as  trial  wave  function,  we 
find  that  the  phaseless  AFQMC  method  generally  gives  better  overall  accuracy  and  a  more 
uniform  behavior  than  the  coupled  cluster  CCSD(T)  method  in  mapping  the  potential-energy 
curve  [12,  13,  14].  Calculations  were  also  done  in  the  triple-bond  system  of  N2,  and  a  chain 
of  50  H  atoms. 


D  Eliminating  spin  contamination  in  the  treatment  of 
magnetic  properties. 

The  use  of  an  approximate  reference  state  wave  function  \^t)  hr  electronic  many-body 
methods  can  break  the  spin  symmetry  of  Born-Oppenheimer  spin-independent  Hamiltonians. 
This  can  result  in  significant  errors.  To  address  the  problem,  a  simple  spin-projection  method 
was  introduced  [13]  for  AFQMC,  which  yields  spin-contamination-free  results,  even  when  a 
spin-contaminated  |Tt)  is  used.  The  method  was  applied  to  the  difficult  F2  molecule,  which 
is  unbound  within  UHF.  With  a  UHF  | T^),  spin  contamination  causes  large  systematic 
errors  and  long  equilibration  times  in  AFQMC  in  the  intermediate,  bond-breaking  region. 
The  spin-projection  method  eliminated  these  problems,  and  delivered  an  accurate  potential 
energy  curve  from  equilibrium  to  the  dissociation  limit  using  the  UHF  |4W).  The  calculated 
spectroscopic  constants,  using  large  basis  sets,  were  in  excellent  agreement  with  experiment. 
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Figure  2:  Spin-density  wave  (SDW)  states  in  the  homogeneous  electron  gas,  in  Hartree-Fock  (HF) 
theory  [19].  The  left  shows  the  result  from  our  full  unrestricted  HF  calculation.  The  HF  ground  state 
clearly  exhibits  SDW  structure,  shown  by  the  spin  density  for  a  plane  in  the  supercell.  The  right 
shows  a  cartoon  depicting  our  theory  explaining  the  mechanism  for  the  SDW  state,  which  we  showed 
to  have  characteristic  wave  vector  qa  <  2 kp,  not  2 kp  as  previously  believed  [20,  21]. 


E  Finite-size  correction  in  many-body  electronic  struc¬ 
ture  calculations 

Finite-size  effects  are  a  major  source  of  error  in  many-body  electronic  structure  calculations 
of  extended  systems.  We  have  developed  a  general  method  to  correct  for  such  errors  [7].  We 
showed  that  many-body  finite-size  effects  can  be  effectively  included  in  a  modified  local  den¬ 
sity  approximation  calculation.  A  parameterization  for  the  finite-size  exchange-correlation 
functional  was  obtained.  The  method  is  simple  and  gives  post-processing  corrections  that  can 
be  applied  to  any  many-body  calculations.  Applications  to  a  model  insulator,  to  semicon¬ 
ducting  Si,  and  to  metallic  Na  showed  that  the  method  delivers  greatly  improved  finite-size 
corrections.  The  method  has  recently  been  implemented  by  Dario  Alfe  in  the  open-source 
PWSCF  package  [17]  (and  used  in  calculations  of  solid  iron  [18].) 

F  What  is  the  Hartree-Fock  ground  state  in  the  elec¬ 
tron  gas? 

In  1962,  Overhauser  showed  that  within  HF  the  electron  gas  is  unstable  to  a  spin  density 
wave  (SDW)  state.  Determining  the  true  HF  ground  state  has  remained  a  challenge.  Using 
numerical  and  analytic  techniques,  we  obtained  the  HF  ground  state  of  the  three-dimensional 
electron  gas.  At  high  density,  we  found  a  number  of  novel  states  characterized  by  nearly 
constant  charge  density  but  broken  spin  symmetry,  with  an  SDW  wave  vectors  differing  from 
the  Overhauser  states.  We  showed  that  the  broken- symmetry  state  originates  from  pairing 
instabilities  at  the  Fermi  surface,  and  proposed  a  model  to  explain  its  properties  [19]  (Fig.  2). 
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Figure  4:  Accurate  calculations  of  the  equation  of  state  (EOS)  in  the  2-D  repulsive  Hubbard 
model.  This  was  made  possible  by  recent  advances  in  using  twist-averaging  boundary  conditions 
(TABC)  [22]  and  in  controlling  the  phase  problem[23,  24,  1],  which  eliminated  shell  effects  (see 
kinks  on  the  left  from  a  typical  prior  calculation  [25])  and  drastically  reduced  finite-size  effects. 
EOS  curves  on  the  right  are  smooth  and  different  lattice  sizes  are  indistinguishable.  A  fit  for  the 
EOS  was  obtained.  From  the  inset  on  the  right,  it  is  seen  that  the  quality  of  the  fit  is  excellent. 


G  Two-dimensional  optical  lattices  with  repulsive  in¬ 
teraction  —  equation  of  state. 

We  have  determined  the  equation  of  state  in  the  ground  state  of  the  two-dimensional  repul¬ 
sive  single-band  Hubbard  model  at  intermediate  interactions.  Shell  effects  are  eliminated 
and  finite-size  effects  are  greatly  reduced  by  boundary  condition  integration.  Spin-spin  cor¬ 
relation  functions  and  structure  factors  are  also  calculated.  In  lattice  sizes  up  to  16  x  16,  the 
results  show  signal  for  phase  instability.  Upon  doping,  the  energy  in  these  finite  systems  is 
larger  than  it  would  be  if  the  system  were  to  separate  into  two  phases:  one  of  density  n  —  1 
(hole- free)  and  the  other  at  density  nc  (~  0.9).  The  long-range  antiferromagnetic  order  is 
coupled  to  this  process,  and  is  lost  below  nc.  This  suggests  that  in  the  thermodynamic  limit, 
either  the  system  truely  phase  separates,  or  there  exists  long  wavelength  collective  modes 
whose  wavelength  is  comparable  to  or  larger  than  the  system  size  we  are  able  to  simulate. 
The  question  is  addressed  in  the  next  part. 


H  Spin-density  wave  states  with  long  wave  length 
modulation  in  the  ground  state  of  the  doped  Hub¬ 
bard  model. 

What  happens  to  the  antiferromagnetic  (AF)  order  in  the  Hubbard  model  when  it  is  doped? 
This  is  a  long-standing  fundamental  question  not  only  for  understanding  the  magnetic  prop¬ 
erties,  which  the  Hubbard  model  was  originally  designed  to  describe,  but  also  in  the  context 
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Figure  3:  SDW  states  with  long  wavelength  modulation  in  the  2-D  Hubbard  model  [26].  The  left 
panel  is  a  3-D  plot  of  the  spin-spin  correlation  function  in  the  ground  state  of  8  x  32  lattice  with 
240  electrons  (n  =  0.9375)  and  U/t  —  4.  Anti-ferromagnetic  order  can  be  seen  from  the  diagonal 
“lines."  Phase  shifts  occur  at  the  nodal  lines,  across  which  there  is  a  mismatch  of  the  diagonal  lines. 
The  right  panel  examines  the  doping  dependence.  The  staggered  spin-spin  correlation  function  is 
plotted,  at  three  different  densities.  Calculations  are  done  using  free-electron  trial  wave  function.  The 
long  wavelength  modulation  and  nodal  lines  (two  colors  erasing)  are  clear.  As  doping  is  increased,  the 
wavelength  of  the  modulating  wave  decreases,  as  does  the  amplitude  of  the  SDW. 


of  high-Tc,  which  shares  the  same  parameter  regime  and  is  believed  to  be  closely  related  to 
AF  fluctuations.  We  have  calculated  the  ground-state  equation  of  state  [24,  26]  and  the  mag¬ 
netic  correlations  [26] ,  using  the  phaseless  AFQMC  method  (which  reduces  to  the  constrained 
path  method  [27]  with  on-site  interaction).  Twist-average  boundary  conditions  were  used, 
a  common  algorithmic  development  with  electronic  structure  calculations  [7].  An  accurate 
equation  of  state  was  obtained  for  the  paramagnetic  phase  (doping  >~  0.15).  With  smaller 
doping,  an  incommensurate  SDW  state  is  found  at  intermediate  interaction  strengths,  with 
antiferromagnetic  order  and  essentially  homogeneous  charge  correlation.  The  wavelength 
of  the  collective  mode  decreases  with  doping,  as  does  its  magnitude.  As  the  interaction  is 
increased,  the  holes  go  from  a  wave-  to  a  particle-like  state,  and  charge  ordering  begins  to 
develop  which  eventually  evolves  into  stripe  states  [28]. 


References 

[1]  Wirawan  Purwanto  and  Shiwei  Zhang.  Correlation  effects  in  the  ground  state  of  trapped 
atomic  bose  gases.  Phys.  Rev.  A,  72(5):053610,  lONovember  2005. 

[2]  Wirawan  Purwanto  and  Shiwei  Zhang.  Quantum  Monte  Carlo  method  for  the  ground  state  of 
many-boson  systems.  Phys.  Rev.  F,  70:056702,  2004. 

[3]  W.  A.  Al-Saidi,  Shiwei  Zhang,  and  Henry  Krakauer.  Auxiliary-field  quantum  Monte  Carlo 
calculations  of  molecular  systems  with  a  Gaussian  basis.  The  Journal  of  Chemical  Physics , 
124(22)  :224101,  2006. 


5 


[4]  Shiwei  Zhang,  Henry  Krakauer,  Wissam  A.  Al-Saidi,  and  Malliga  Suewattana.  Quantum 
simulations  of  realistic  systems  by  auxiliary  fields.  Comput.  Phys.  Commun .,  169:394,  2005. 

[5]  W.  A.  Al-Saidi,  Henry  Krakauer,  and  Shiwei  Zhang.  Auxiliary-field  quantum  Monte  Carlo 
study  of  TiO  and  MnO  molecules.  Physical  Review  B  ( Condensed  Matter  and  Materials 
Physics ),  73(7):075103,  2006. 

[6]  Shiwei  Zhang  and  Henry  Krakauer.  Quantum  Monte  Carlo  method  using  phase-free  random 
walks  with  Slater  determinants.  Phys.  Rev.  Lett.,  90(13) :  136401 ,  Apr  2003. 

[7]  Hendra  Kwee,  Shiwei  Zhang,  and  Henry  Krakauer.  Finite-size  correction  in  many-body  elec¬ 
tronic  structure  calculations.  Physical  Review  Letters ,  100(12):126404,  2008. 

[8]  W.  A.  Al-Saidi,  Henry  Krakauer,  and  Shiwei  Zhang.  Auxiliary-field  quantum  Monte 
Carlo  study  of  first-  and  second-row  post-d  elements.  The  Journal  of  Chemical  Physics , 
125(15) :  1541 10,  2006. 

[9]  W.  A.  Al-Saidi,  Henry  Krakauer,  and  Shiwei  Zhang.  A  study  of  H  +  H2  and  several  H-bonded 
molecules  by  phaseless  auxiliary-field  quantum  Monte  Carlo  with  plane  wave  and  Gaussian 
basis  sets.  The  Journal  of  Chemical  Physics,  126(19):194105,  2007. 

[10]  Rodney  J.  Bartlett  and  Monika  Musial.  Coupled-cluster  theory  in  quantum  chemistry.  Rev. 
Mod.  Phys.,  79(1):291,  2007. 

[11]  T.  Daniel  Crawford  and  Henry  F.  Schaefer  III.  An  introduction  to  coupled  cluster  theory  for 
computational  chemists.  Reviews  in  Computational  Chemistry,  14:33-136,  2000. 

[12]  W.  A.  Al-Saidi,  Shiwei  Zhang,  and  Henry  Krakauer.  Bond  breaking  with  auxiliary-field  quan¬ 
tum  Monte  Carlo.  The  Journal  of  Chemical  Physics,  127(14):144101,  2007. 

[13]  Wirawan  Purwanto,  W.  A.  Al-Saidi,  Henry  Krakauer,  and  Shiwei  Zhang.  Eliminating  spin 
contamination  in  auxiliary- field  quantum  Monte  Carlo:  Realistic  potential  energy  curve  of  F2. 
The  Journal  of  Chemical  Physics,  128(11):  114309,  2008. 

[14]  Wirawan  Purwanto,  Shiwei  Zhang,  and  Henry  Krakauer.  Excited  state  calculations  using 
phaseless  auxiliary- field  quantum  Monte  Carlo:  Potential  energy  curves  of  low-lying  C2  singlet 
states.  The  Journal  of  Chemical  Physics,  130(9) :094107,  2009. 

[15]  Malliga  Suewattana,  Wirawan  Purwanto,  Shiwei  Zhang,  Henry  Krakauer,  and  Eric  J.  Walter. 
Phaseless  auxiliary-field  quantum  Monte  Carlo  calculations  with  plane  waves  and  pseudopo¬ 
tentials:  Applications  to  atoms  and  molecules.  Physical  Review  B  (Condensed  Matter  and 
Materials  Physics),  75(24) :245123,  2007. 

[16]  K  P  Esler,  J  Kim,  D  M  Ceperley,  W  Purwanto,  E  J  Walter,  H  Krakauer,  S  Zhang,  PRC  Kent, 
R  G  Hennig,  C  Umrigar,  M  Bajdich,  J  Kolorenc,  L  Mitas,  and  A  Srinivasan.  Quantum  Monte 
Carlo  algorithms  for  electronic  structure  at  the  petascale;  the  Endstation  project.  Journal  of 
Physics:  Conference  Series,  125:012057  (15pp),  2008. 

[17]  S.  Baroni,  A.  Dal  Corso,  S.  de  Gironcoli,  and  P.  Giannozzi.  http://www.pwscf.org. 

[18]  Ester  Sola  and  Dario  Alfe.  Melting  of  iron  under  earth’s  core  conditions  from  diffusion  monte 
carlo  free  energy  calculations.  Physical  Review  Letters,  103:078501,  2009. 


6 


[19]  Shiwei  Zhang  and  D.  M.  Ceperley.  Hartree-Fock  Ground  State  of  the  Three-Dimensional 
Electron  Gas.  Phys.  Rev.  Lett.,  100:236404,  2008. 

[20]  A.  W.  Overhauser.  Spin  density  waves  in  an  electron  gas.  Phys.  Rev.,  128:1427,  1962. 

[21]  G.  F.  Giuliani  and  G.  Vignale.  Quantum  Theory  of  the  Electon  Liquid.  Cambridge  University 
Press,  Cambridge,  UK,  2005. 

[22]  C.  Lin,  F.-H.  Zong,  and  D.  M.  Ceperley.  Twist-averaged  boundary  conditions  in  continuum 
Quantum  Monte  Carlo  algorithms.  Phys.  Rev.  E ,  64:016702,  2001. 

[23]  Shiwei  Zhang  and  Henry  Krakauer.  Quantum  Monte  Carlo  method  using  phase-free  random 
walks  with  Slater  determinants.  Phys.  Rev.  Lett.,  90:136401,  2003. 

[24]  Chia-Chen  Chang  and  Shiwei  Zhang.  Spatially  inhomogeneous  phase  in  the  two-dimensional 
repulsive  Hubbard  model.  Physical  Review  B ,  78:165101,  2008. 

[25]  Nobuo  Furukawa  and  Masatoshi  Imada.  Two-Dimensional  Hubbard  Model  Metal  Insulator 
Transition  Studied  by  Monte  Carlo  Calculation.  J.  Phys.  Soc.  Jpn .,  61,  1992. 

[26]  Chia-Chen  Chang  and  Shiwei  Zhang.  Spin  and  charge  order  in  doped  Hub¬ 
bard  model:  long- wavelength  collective  modes.  Phys.  Rev.  Lett,  (submitted),  2009. 

http:/ /arxiv.org/abs/0910.0712. 

[27]  Shiwei  Zhang,  J.  Carlson,  and  J.  E.  Gubernatis.  Constrained  Path  Monte  Carlo  Method  for 
Fermion  Ground  States.  Phys.  Rev.  B,  55:7464,  1997. 

[28]  S.  R.  White  and  D.  J.  Scalapino.  Stripes  on  a  6-leg  hubbard  ladder.  Phys.  Rev.  Lett.,  91:136403, 
2003. 


7 


